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ABSTRACT: The motivation, present status, and future plans of the search for f3/3(0u) decay are reviewed. It is 
argued that, motivated by the recent observations of neutrino oscillations, there is a reasonable hope that f5/3(0v) decay 
corresponding to the neutrino mass scale suggested by oscillations, m„ ~ 50 meV, actually exists. The challenges to 
achieve the sensitivity corresponding to this mass scale, and plans to overcome them, are described. 
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1 INTRODUCTION 

Since the last Annual Review article on double beta decay (||), published in 1994, there have been 
several exciting developments. Most significantly, the neutrino oscillation experiments convincingly 
show that neutrinos have a finite mass. However, in oscillation experiments only the differences in 
squares of the neutrino masses, Am| = |m| — mj\, can be measured. Nevertheless, a lower limit 
on the absolute value of the neutrino mass scale, m sca i e = V Am 2 , has been established in this way. 
Its existence, in turn, is causing a renaissance of enthusiasm in the double beta decay community 
which is expected to reach, in the next generation of experiments, the sensitivity corresponding to 
this mass scale. Below we review the current status of the double beta decay and the effort devoted 
to reach the required sensitivity. But before proceeding, we briefly summarize the achievements of 
the neutrino oscillation searches and the role that the search for the neutrinoless double beta decay 
plays in the elucidation of the pattern of neutrino masses and mixing. In these introductory remarks 
we use the established terminology, some of which will only be defined later in the text. 

There is a consensus that the measurement of atmospheric neutrinos by the SuperKamiokande 
collaboration (|2|) can be only interpreted as a consequence of the nearly maximum mixing between 
and v T neutrinos (i.e., sin 2 29 atm ~ 1.0) with the corresponding mass squared difference Am 2 tTO ~ 
3 x 10~ 3 eV 2 , i.e., m sca i e ~ 50 meV. This finding is supported by the K2K experiment (|3|) that uses 
an accelerator v„ beam pointing towards the SuperKamiokande detector 250 km away, and is in 
accord with the earlier findings of the Kamiokande (|j), 1MB (||) and Soudan (|6|) experiments. Several 
large long-baseline experiments are being built to further elucidate this discovery, and determine the 
corresponding parameters more accurately. 

At the same time the 'solar neutrino puzzle', which has been with us for over thirty years since the 
pioneering chlorine experiment of Davis (0), also reached the stage where the interpretation of the 
measurements in terms of oscillations between the v e and some combination of active, i.e. and v T 
neutrinos, seems inescapable. In particular, the juxtaposition of the results of the SNO experiment 
(||) and SuperKamiokande (^), together with the earlier solar neutrino flux determination in the 
chlorine and gallium (1^,11) experiments, leads to that conclusion. The value of the corresponding 



oscillation parameters remain uncertain, with several 'solutions' possible, although the so-called Large 
Mixing Angle (LMA) solution with sin 2 2Q soi ~ 0.8 (but sin 2 26 soi < 1) and Am 2 sol ~5x 10~ 5 eV 2 is 
preferred at present. Again, the continuing and soon to be operational experiments, like KamLAND 
and Borexino, aim to find with more certainty which of the possible solutions is the correct one. 
The pattern of neutrino mixing is further simplified by the constraint due to the Chooz and Palo 



Verde reactor neutrino experiments (12, 13) which lead to the conclusion that the third mixing angle, 
013, is small, sin 2 2#i 3 < 0.1. 

The oscillation experiments cannot determine the absolute magnitude of the masses and, in partic- 
ular, cannot at this stage separate two rather different scenarios, the hierarchical patt ern of neutrino 
masses in which the neutrino masses mj and/or rrij are of similar magnitude as J Am| and the 

degenerate pattern in which all rrii S> ^jAmfj . It is likely that the search for the neutrinoless double 
beta decay, reviewed here, will help in the foreseeable future in establishing the correct mass pattern 
and in determining, or at least strongly constraining, the absolute neutrino mass scale. 

Moreover, the oscillation results do not tell us anything about the properties of neutrinos under 
charge conjugation. While the charged leptons are Dirac particles, distinct from their antiparticles, 
neutrinos may be the ultimate neutral particles, as envisioned by Majorana, identical with their an- 
tiparticles. That fundamental distinction becomes important only for massive particles. Neutrinoless 
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double beta decay proceeds only when neutrinos are massive Majorana particles, hence its observation 
would resolve the question. 

Quite generally, the Standard Electroweak Model postulates that neutrinos are massless, and that 
the total lepton number, as well as the individual flavor lepton numbers, are conserved. For various 
reasons most people believe that the model, despite its enormous success, is incomplete, and the hunt 
for 'physics beyond the Standard Model' is being actively pursued on many fronts. Observation of 
neutrino mass and oscillation is clearly an example of a phenomenon that is at variance with the 
Standard Model. Further elucidation of various aspects of neutrino mass will undoubtedly help in 
pointing to the proper generalization of the Standard Model. 

Double beta decay is a rare transition between two nuclei with the same mass number A involving 
change of the nuclear charge Z by two units. The decay can proceed only if the initial nucleus is 
less bound than the final one, and both must be more bound than the intermediate nucleus (or the 
decay to the intermediate nucleus must be highly hindered, as in 48 Ca). These conditions are fulfilled 
in nature for many even-even nuclei, and only for them. Typically, the decay can proceed from the 
ground state (spin and parity always + ) of the initial nucleus to the ground state (also + ) of the 
final nucleus, although the decay into excited states (0 + or 2+) is in some cases also energetically 
possible. 

The two-neutrino decay, (5(5(2v), 

(Z,A) -> (Z + 2,A) + ef + eg + v e i + v e2 (1) 

conserves not only electric charge but also lepton number. (Analogous decays, involving transforma- 
tion of two protons into two neutrons, are also sometimes possible. We concentrate here, however, 
on the decays 2n — > 2p with more candidate nuclei, and usually larger Q values.) 
On the other hand, the neutrinoless decay, (5(5(0i/), 

(Z,A) -» (Z + 2,A) + (2) 

violates lepton number conservation and is therefore forbidden in the standard electroweak theory. In 
addition, there can be transitions (5(5(0v, x) m which a light neutral boson x, a Majoron postulated 
in various extensions of the standard electroweak theory fll4|), is emitted: 

(Z,A) -^{Z + 2,A) + e^ + e^+ X • (3) 

It is fascinating to realize that the interest in (5(5 decay spans more than six decades. Already in 
1937 Racah (|l^), following the fundamental suggestion of Majorana fll6|), discussed the possibility 
of a neutrinoless transformation of two neutrons into two protons plus two electrons. Even earlier, 
Goeppert-Mayer (|l7|) evaluated the decay rate of the (5(5{2v) mode and realized that the corresponding 
half-life could exceed 10 20 y. And Furry, shortly afterwards fll8|) , estimated that the (3/3(0v) decay 
should be much faster than the (5(5{2v) decay. That conclusion, however, had to be revised with the 
discovery of parity nonconservation in weak interactions. Thus, the stage was set for the realization 
that the observation of the (5(5(0u) decay would establish that the neutrino is a massive Majorana 
particle. 

There were numerous earlier reviews of (5(5 decay, including the classics by Primakoff & Rosen 
(|l9|), Haxton & Stephenson (f2C|), and Doi, Kotani & Takasugi (|2l|). More recent reviews, besides 
([[]), include those by Boehm & Vogel ( |22| ) (where the phase-space integrals are listed), Suhonen & 
Civitarese (^) (which deals mostly with the nuclear matrix elements), Faessler & Simkovic (pi|) , 
Vergados (p5|), and Klapdor-Kleingrothaus (^). A rather complete list of experimental data can be 
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found in the papers by Tretyak and Zdesenko (27, 28). In the Review of Particle Physics (^) the 
most important (5(5 decay experimental data are regularly listed. The whole field of neutrino mass 
and oscillations has been reviewed recently in Annual Reviews (|30"1). 

It is easy to distinguish the three decay modes ((5(5(2v), /9/3(Oi/), and (5(5(0v, x)) by the shape 
of the electron sum energy spectra, which are determined by the phase space of the outgoing light 
particles. For the (5(5{2v) and (5(3(0i>) modes these spectra are illustrated in Figure [l]. In the 2v decay 
there is a broad maximum at the sum kinetic energy of the two electrons below half of the endpoint 
energy. In contrast, in the (5(3(0is) mode the two electrons carry the full available kinetic energy (the 
nuclear recoil is negligible for all practical purposes) and the spectrum is therefore a single peak at the 
endpoint. In the Majoron decay, not shown in order not to clutter the figure, the electron spectrum 
is again continuous, but the maximum is shifted higher, above the halfway point, as required by the 
three-body light particle phase space. 

The insert in Fig. [l] illustrates in detail the expected spectra near the endpoint where the (5(5(2v) 
decay represents the ultimate background in the search of the /?/5(0z^) mode. (See section 4.2.4 for 
the discussion of this point.) 

The (5(5{2v) decay mode is an allowed process. However, since it is a second order semileptonic weak 
decay, its lifetime, proportional to (Gf cos 0c)~ A , is very long. The observation of the (5(5{2v) decay 
presents a formidable challenge, since it must be detected despite the presence of inevitable traces 
of radioisotopes with similar decay energy, but lifetimes more than 10 orders of magnitude shorter. 
Nevertheless, at present, that challenge has been met and the (5(5{2v) decay has been positively 
identified in a number of cases. Observing the (5(5{2i/) decay is important not only as a proof that 
the necessary background suppression has been achieved, but also allows one to constrain the nuclear 
models needed to evaluate the corresponding nuclear matrix elements. 

Hand in hand with the observation of the (5(5{2v) decays, the experiments became sensitive to 
longer and longer half-lives of the (5(5{Qv) decay mode. Since the rate of (5(5{Qu) is proportional to 
the square of effective neutrino mass, the improvements lead to correspondingly improved limits on 
the mass. This is illustrated in Figure ^ which shows an essentially logarithmic improvement, by 
a factor larger than about four every decade, of the corresponding limits. If this trend continues, 
we expect to reach the neutrino mass scale suggested by the oscillation experiments in 10-20 years. 
Given the typical lead time of the large particle physics experiments, it suggests that the relevant (5(5 
decay experiments should enter the 'incubation' process now. 



2 NEUTRINO MASS - THEORETICAL ASPECTS 



2.1 Majorana and Dirac neutrinos 

Empirically, neutrino masses are much smaller than the masses of the charged leptons with which 
they form weak isodoublets. Even the mass of the lightest charged lepton, the electron, is at least 
10 5 times larger than the neutrino mass constrained by the tritium beta decay experiments. The 
existence of such large factors is difficult to explain, unless one invokes some symmetry principle. 
The assumption that neutrinos are Majorana particles is often used in this context. Moreover, many 
theoretical constructs invoked to explain neutrino masses lead to this conclusion. 

The term 'Majorana' is used for particles that are identical with their own antiparticles while Dirac 
particles can be distinguished from their antiparticles. This implies that Majorana fermions are two- 
component objects, while Dirac fermions are four-component. In order to avoid confusion and to 
derive the formula for the /3/?(0z^) rate mediated by the exchange of massive Majorana neutrinos, it 
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is worthwhile to discuss, albeit briefly, the formalism needed to describe them (for more details, see 



e.g. (|31|,|21|,|32|,|3 



Massive fermions are usually described by the Dirac equation, where the chirality eigenstates ipR 
and if)L are coupled and form a four-component object of mass m, 

i^d^ijjR - rmp L = , i(a^d^ L - m^ R = , (4) 

where = (a , a), = (a ,— a) and (cr°,a) are the Pauli matrices. As written, ipi(R) are two- 
component spinors; the usual four-component bispinors are defined as: 

where ^l(r) are just the chiral projections of i.e. the eigenstates of Pur) = (1 T 7s)/2- 



However, Majorana's suggestion (16) allows one to use an alternative description of those massive 
fermions which do not have any additive quantum numbers as either two-component ipR (mass m), 
or ipL (mass m'), which obey independent equations 

Wfl - me^R = ; i{^3^ L + m'e^l = , (6) 

where e = ia y . 

The Majorana fields can be also expressed in the four-component form 

= ( ~t t ' ) ' * rW = ( ) ' (7) 

Such a four-component notation is a convention useful to express the charged weak current in a 
compact form. It is then clear that the Dirac field Vl/, Equation ||, is equivalent to a pair of Majorana 
fields with m = m' and ipL = eip R . 

The four-component Majorana fields, Equation [7|, are selfconjugate, = *&l(r)( x )i where 

charge conjugation is defined as ^ L m\(x) = ^ 2 ^°^'l(R)- ^ ne fields ^l{x) and ^r(x) are eigenstates 
of CP with opposite eigenvalues. 

The Lorentz invariant mass term in the neutrino Langrangian can appear in three forms: 

M d [9rv l + (vl) c v c r\ , M L \{p£fv L + v L v c l] , Mr[{vr) c ur + uru c r } , (8) 

where we have introduced the notation vl(R) f° r t ne corresponding neutrino annihilation operators. 
The first expression in Equation |8| is the Dirac mass term (with the mass parameter Md) which 
requires the existence of both chirality eigenstates vl and vr and conserves the lepton quantum 
number. The second (and third) mass terms are Majorana mass terms, which violate the lepton 
number and can be present even without the existence of vr (for the term with mass parameter Ml) 
or vl (for the term with mass parameter Mr). In general, all three terms might coexist, and then the 
mass Langrangian must be diagonalized resulting in two generally nondegenerate mass eigenvalues 
for each flavor. (That is the situation with the generic see-saw mass (|34j), where it is assumed that 
Mr 3> Md 3> Ml ~ 0, and the light neutrino acquires the mass m v ~ M^/Mr.) 

Let us consider now the general situation with N flavors of the left-handed neutrinos vl and in 
addition an equal number ./V of the right handed neutrinos vr. The most general Lorentz invariant 
mass term of the neutrino Langrangian has then the form 
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where vl and vr are column vectors of dimension N. Here Ml and Mr are symmetric N x N 
matrices (Majorana masses for the left- and right-handed neutrinos) and Mb is an arbitrary and 
generally complex N x N matrix. 

The mass matrix M, with real positive eigenvalues mi, . . . ,m2Ni is diagonalized by the 2N x 2N 
unitary matrix 




The general mixing matrices U and V have N rows and 2N columns and is a column vector of 
dimension 2N of Majorana-like objects (pip. On the other hand, if none of the states Ur exist, or 
if Mr is so large that the corresponding states need not be considered, only Ml is relevant, and 
only the N x N mixing matrix U is needed to diagonalize the mass term (and <&l nas then only N 
components, naturally). 

Let us consider in more detail the latter case, when only Ml 7^ 0. The N x N unitary mixing 
matrix U contains N 2 real parameters. However, N of them correspond to unphysical phases; there 
are N(N — l)/2 angles and N(N — l)/2 physically relevant phases describing possible CP violations. 
(For a discussion of parameter counting, see (pl|) .) In the oscillation experiments that violate only 
the flavor lepton number, but conserve the total lepton number (such as u e — > ia, or i/„ — > u T ), one 
can determine, in principle, all angles and (N — 1)(N — 2)/2 phases. These phases, common to the 
Dirac and Majorana neutrinos, describe CP violation responsible for the possible differences of the 
oscillation probabilities i>t — > vy and &t — * . 

The remaining N — 1 phases affect only neutrino oscillation-like processes (in which neutrinos are 
created in the charged current weak processes and absorbed again in charged current) that violate 
the total lepton number, such as the (3/3 (0v) decay. Such phases are physically significant only for 
Majorana neutrinos; they are unphysical for Dirac neutrinos. This is so because for Majorana neutri- 
nos one cannot perform the transformation V{ — > u[ = e lOLi Vi, which would violate the selfconjugation 
property. 

In principle the distinction between Dirac and Majorana neutrinos affects other processes as well, 
such as the angular distribution of v — e scattering or photon polarization in the — > Uj + 7 decay. 
However, the "Practical Dirac-Majorana Confusion Theorem" ( |35| ) states that the distinction vanishes 
for m u — > which makes it essentially unobservable in these cases. 



2.2 PP(0u) decay rate and Majorana mass 

Here we shall consider only the simplest case of the left-handed V — A weak currents and light massive 
Majorana neutrinos. This is the case of current interest provided the neutrino mass revealed in the 
oscillation experiments is of Majorana character. The more general expressions can be found e.g. in 
the reviews (^,^Tj). (For recent formulation of the general problem, see fl36|).) 
The differential decay rate of the (3(3(Qv) process, Equation g is @ 



aT 0u = 2nJ2 l«o,| 2 5(ei + e 2 + E f - M,)^ ^ , (11) 



spin 



where 61(2) and f>\(2) are total energies and momenta of the electrons and Ef{Mi) is the energy of the 
final (mass of the initial) nuclear state. The quantity Rq u is the reaction amplitude to be evaluated 
in the second order pertubation theory with respect to the weak interactions. 
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The lepton part of Rov, involving the emission and reabsorption of the Majorana neutrino of mass 
nij , is 

where Pl = (1 — 7s)/2, e(x), e c (y) are the electron creation operators, and q is the momentum transfer 
four- vector. Since 7^ anticommute with 75, this amplitude is proportional to rrij and the term with 
9^7/i vanishes. After integrating over the energy of the virtual neutrino dq°, the denominator q 2 — m 2 

is replaced by its residue uij/ir, where uij = y 1 q 2 + m 2 - . The amplitude is therefore proportional to 
rrij/ujj <C 1 for light neutrinos. 

The remaining integration over the virtual neutrino momentum q leads to the appearance of the 
neutrino potentials 

H k (r, A k ) = ™*L f°° dq ^^J, , A m = E m - (M t + M f )/2 ± (ei - e 2 )/2 , (13) 
irr Jo uj(lo + Ak) 

where 1 and 2 label the emitted electrons, E m is the excitation energy of the intermediate nucleus, 
Mf is the mass of the final nucleus, and r is the distance between the two neutrons that are changed 
into protons. The factor Rn, the nuclear radius, is introduced in order to make the potential H 
dimensionless. In the case of the (5(3{Qv) decay one can use the closure approximation, replacing E m 
by an appropriate mean value. (This is justified because we expect that the momentum of the virtual 
neutrino is determined by the uncertainty relation q ~ 1/r ~ 100 MeV, thus the variation of E m from 
state to state can be neglected.) The contributions of the two electrons are then added coherently, 
and thus the neutrino potential to use is 

H(r) = [H 1 (r,A 1 ) + H 2 (r,A 2 )]/2^H(r,A) , (14) 

where A = E m — (Mi + M/)/2 and E m is the average energy of the intermediate nucleus. The 
potential H(r) only very weakly depends on rrij as long as the neutrino mass is less than ~ 10 MeV. 

For the ground state to ground state, i.e., 0^ — > 0j transitions, it is enough to consider s-wave 
outgoing electrons, and the nonrelativistic approximation for the nucleons. The nuclear part of the 
amplitude then turns into a sum of the Gamow- Teller and Fermi nuclear matrix elements, where the 
superscript Of is used to signify the presence of the neutrino potential H (r): 

\M 0v \ = M°J T - 4 M F = (f\ E H (nk, A)t+t+ U -S k - 9 \\ \i) . (15) 
9a it \ 9aJ 

The summation is over all nucleons, |/) are the final (initial) nuclear states, and gv(dA) are 
the vector (axial vector) coupling constants. Such an expression is now analogous to the allowed 
approximation of the ordinary beta decay. 

Thus, in the approximations described above, which are quite accurate, the transition amplitude 
for a Majorana neutrino of mass rrij is simply a product of rrij and the above combination of the 
nuclear matrix elements. However, since in each of the two vertices an electron is emitted, the mixing 
amplitude U e j apears in each of them, and the physical (5f3(Qu) reaction amplitude contains the factor 
U 2 ,j (not |t/ e j| 2 ) and is proportional to the factor 

(16) 



{m v ) 
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where the sum is only over light neutrinos with rrij < 10 MeV (for heavier ones one cannot neglect 
the mass in the neutrino propagator, Equation |l^). The quantity (m„) is the effective neutrino mass. 
Since and not \U e j\ 2 appear in (m u ), its value depends on the Majorana phases discussed above. 
To obtain the decay rate, the reaction amplitude has to be squared, and multiplied by the corre- 



sponding phase space integral, which in this case, see Equation 11, is simply the two-electron phase 
space integral proportional to 



G 0u ~ J F(Z,e 1 )F(Z,e 2 ) P iP2eie 2 5(Eo - e x - e 2 )de l de 2 , 



(17) 



where Eq is the available energy (the sum electron kinetic energy peak is at Q = Eq — 2m e ). F(Z, e) 
is the usual Fermi function that describes the Coulomb effect on the outgoing electron. 

Summarizing, if the /3/3(Oz/) decay is mediated by the exchange of a light massive Majorana neutrino 
(the assumption that we wish to test), the half-life is 

2 



H 



G°»(E 0: Z) 



M°Jt 



& M 0u 
2 Nip 

9a 



(18) 



where G 0u is the exactly calculable phase space integral, (m u ) is the effective neutrino mass and Mq Tj 
M F V are the nuclear matrix elements, defined in Equation [D| The way these nuclear matrix elements 
are evaluated, and the associated uncertainty, is discussed in the next Section. (As explained earlier, 
the neutrino mass appears in the amplitude in the combination nij/ujj <C 1; the denominator uij has 
been absorbed in the neutrino potential H(r).) 

Thus, if an upper limit on f3(3(0v) rate is experimentally established, and the nuclear matrix 
elements are known, one can deduce the corresponding upper limit on {m u ). On the other hand, 
if (3(3(Qv) is observed, one can deduce the appropriate value of (m u ). That is a justified procedure, 
however, only if the exchange of the light Majorana neutrino, discussed above, is indeed the mechanism 
responsible for the decay. There in no way to decide on the mechanism when only the decay rate 



is known. However, a general theorem (37) states that once (3(3(Qv) has been observed, in gauge 
theories the Majorana neutrino mass necessarily arises. But the magnitude of the corresponding 
neutrino mass is difficult to estimate if the exchange of a virtual light Majorana neutrino is not the 
dominant mechanism of the f3f3(0u) decay. 



2.3 /3/3(Oz/) decay and oscillation parameters 

Let us assume that there are N massive Majorana neutrinos z^, i = 1, . . . , N. In that case the weak 
eigenstate neutrinos v e ,Vn and v T can be expressed as superpositions of using the 3 x N mixing 
matrix [%. In particular, electron neutrinos are then superpositions, 

N 

"e = Yl U ™ U i ' ( 19 ) 



and the rate of the 0f3(Ov) decay is proportional to (see Equation 18 and Refs. (pq. 39)) 





N 


2 


N 


2 


{mu? = 


U ei m i 

i 




\U ei \ 2 e a 'mi 

i 


, (all mi > 0) 



(20) 



This quantity depends, as indicated, on the N — 1 Majorana phases a.i/2 of the matrix U discussed 
in subsection 2.1 which are irrelevant in neutrino oscillation experiments that do not change the total 
lepton number. 
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If CP is conserved on = kir, but generally any values of 014 are possible. Thus, (m u ) could be 
complex and cancellations in the sum are possible. (For example, a Dirac neutrino corresponds to 
a pair of degenerate Majorana neutrinos with e ai = ±1 whose contribution to {m u ) exactly cancel. 
More generally, some models, such as the Zee model fl40|), postulate that (m u )= 0.) 

While the quantity {m v ) depends on the unknown phases oij, the upper and lower limits of (m„), 

yflv) max 

and (Tn v } m in, depend only on the absolute values of the mixing angles (pl|), 

(m v )max = \ u ei\ 2 fni , {m v ) min = max[(2| U ei \ 2 rrn - {m u ) max ),Q\ . (21) 

i 

Thus, if the search for (3(3(0v) is successful and the value of (m v ) is determined, and at the same time 
the mixing angles \ U e i\ 2 and the mass square differences Am?- are known from oscillation experiments, 
a range of absolute values of the neutrino masses can be deduced. This is illustrated in Figure || 
where we assumed that N = 3, that the Large Mixing Angle (LMA) solution of the solar neutrinos 
is correct, and that the atmospheric neutrino problem requires maximum mixing of the fi and t 
neutrinos. Note that we have to consider two possibilities, the normal and inverted hierarchies (see 
the inserts in Figure |3|) because given the information, we cannot distinguish between them. (Note 
that the uncertainty in the mixing parameters is not included in Figure ||.) 

Naturally, if another constraint exists, for example a successful determination of the neutrino 
mass square J2i Wei\ 2 m 2 in the tritium decay experiments, one can use the knowledge of (m u ) to 
determine or constrain the phases on. 

There have been numerous analyses of the existing data that correlate the current results of the 
neutrino oscillation searches, (3f3(0v) experiments, tritium beta decay experiments, etc. We list here 
only a subset of the corresponding papers, and apologize for omission of other work (see (41, 4^, 43 
H,||,|4|,|3)). 



Altogether, one cannot predict, in general, what the value of (m u ) ought to be using the present 
knowledge. On the other hand, as shown in Figure |3| for the currently most likely oscillation scenario, 
one can show that certain classes of solutions, such as the inverted hierarchy, or the normal hierarchy 



with the smallest neutrino mass 3> y Am 2 ol (degenerate neutrino spectrum) lead to potentially 
observable 0(3(Ov) decay. 

2.4 0/3 (Ou) decay and other lepton number violating processes 

The f3(3{Qv) decay is not the only possible observable manifestation of lepton number violation. 
Muon-positron conversion, 

IT + (A, Z) -> e+ + (A, Z - 2) , (22) 
or rare kaon decays K^^, K een and K^, 

K + -» fi + fi + 7r- , K + -> e + e + vr" , K + -> u + e + ^- , (23) 

are examples of processes that violate total lepton number conservation and where good limits on the 



corresponding branching ratios exist. (See Ref. (48) for a more complete discussion. 



Like the f3/3(Qv) decay, these processes can be mediated by the exchange of a virtual massive 
Majorana neutrino. In that case their rate is proportional to the quantity analogous to (m u ) , 



(m X y) — ^ U x iUyirrii , 

i 



(24) 
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involving other lines of the neutrino mixing matrix U (in (3/3(0v) the relevant quantity is (m ee )). Note 
that, again like in the (3/3(0v), other mechanisms are possible and might lead to faster rates. 

Study of such decays, in principle, would allow one to constrain, or determine, the otherwise 
inaccessible Majorana phases in U. However, the present and foreseeable future of the experimental 
search has not reached the required sensitivity. 

Considerable experimental effort has been devoted to the study of the \i~ — * e + conversion, with 
the best limit © 

T{ ™ + f ~ 6+ + ^ < 1-7 x 10- 12 (90% CL) , (25) 
r(Ti + fi —> + Sc) 

and a substantial improvement is anticipated in proposed experiments. That branching ratio limit 
can be expressed as the limit on 

(m Me ) = (]T U ei U^mi) < 17(82)MeV , (26) 

i 

where the two limiting values reflect the dependence on the spin state (0 or 1) of the created proton 
pair. 

Similarly, for the decay the branching ratio is presently restricted to (see Ref. (^) where the 

limits on the other lepton number violating K decays are described), 

^ < 3.0 x 10- 9 (90% CL) . (27) 



Following Ref. (51), this branching ratio limit can be expressed as the limit 

("*«*) = E U i m i < 4xl0 4 MeV. (28) 

i 

Since, obviously, (m xy ) < m max , the present limits are far from constraining U for the three 
light neutrinos whose mass is restricted by the tritium (3 decay experiments (^,^3)) and the positive 
results of the atmospheric and solar neutrino oscillation experiments to be < O(eV). The hypothetical 
subdominantly coupled heavy sterile are also only marginally constrained (|5l|). 

Nevertheless, it is important to pursue searches for the |AL| = 2 processes with |AL M | ^ 0, since 
they can, generally, yield nonvanishing results even when f3f3(0is) decay is vanishing or very slow. 

3 P/3(0u) MATRIX ELEMENTS 

For all three modes of/?/? decay (/?/?(0f), 0f3(2v), Pf3(0u,x)) one can separate, essentially without loss 
of accuracy, the phase space and nuclear parts of the rate formulae. All nuclear structure effects are 



then represented by the nuclear matrix elements. In the review (23) one can find a rather complete 
list of references, and the results of calculation of the nuclear matrix elements. 

The half-life for the /3/3(2v) decay mode can be written in the compact form, analogous to Equation 



18 but without the factor (m p 



[T 2 ; 2 (0 + - 0+)]- 1 = G 2 »(E ,Z)\M?r T \ 2 , (29) 

where again G 2u is the exactly calculable phase space integral containing all the relevant constants, 
and Mq^ is the nuclear matrix element (there is no Fermi part, due to the isospin conservation), 

MGT - \ E m -(M t + M f )/2 ' (30) 
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where \f)(\i)) are the + ground states of the final (initial) even-even nuclei of masses Mj(Mj), and 
\m) are the 1 + states in the intermediate odd-odd nucleus of energy E m . The individual factors in 
Equation ^ have straightforward physical meaning; the last factor in the numerator is the amplitude 
of the (3~ decay (or of the forward angle (p, n) reaction) of the initial nucleus, while the first factor is 
the amplitude of the (3 + decay (or of the (n,p) reaction) of the final nucleus. Thus the description of 
the (5(5{2v) is equivalent to the description of the full beta strength functions of both the initial and 
final nuclei. The P(3(2v) rate is sensitive to details of nuclear structure, however, because the ground 
state to ground state transition exhausts only a very small fraction of the double GT sum rule (|54j) . 
Description of the (5(3{2v) thus represents a severe test of the nuclear models used in the evaluation 
of M"q T . We list the experimentally determined matrix elements for the (3f3(2v) decay in Table |]. 

10 in the matrix elements, i.e. 



Note that the nuclear structure effects cause variations by a factor 
by a factor ~ 100 in the half-lives. 



The nuclear matrix elements defined in Equation 15 govern both (3f3(Qv) and (3f3(0v, x) decay 
modes. However, the half-life of the (3f3{Qv, \) mode depends on the effective majoron- neutrino 
coupling constant {g vx )-> instead of (m v ) 



0+)] 



-i 



M 



Qu 
GT 



2 M F 
9a 



(fJ, 



v%l i 



(31) 



where again G 0u ' x (Eq, Z) is the phase-space integral, tabulated e.g. in Ref. (21). 

Throughout, we discuss only the PP(0u) decay mediated by the exchange of a light massive Ma- 
jorana neutrino and governed by the nuclear matrix elements (15). The rate of (5(5(Qv) mediated by 
other mechanisms, e.g. involving the right-handed current weak interactions, or the exchange of only 
heavy particles (heavy neutrinos, SUSY particles, etc.) depends on other nuclear matrix elements. 
Detailed discussion of their evaluation is beyond the scope of this review, but an interested reader 
can find it e.g. in Refs. (g|, |||) . 

There are two basic approaches to the evaluation of the nuclear matrix elements for both the (3[3{2u) 
and (3f3{Qv) decays, the quasiparticle random phase approximation (QRPA), and the nuclear shell 
model (NSM). As pointed out above, the factors entering the (3(3{2v) matrix elements are related 
to other nuclear phenomena, and thus testable. This is not so, or at least it is much more difficult, 
for the matrix elements of the /?/3(0z^) mode. It is therefore less clear how to reliably estimate the 
uncertainty involved in their evaluation. 

QRPA has been the most popular theoretical tool in the recent past since it was able to explain the 
suppression, relative to the sum rule (|54|), of M"q T , and it is easy to use. The main physics ingredients 
of the method, relevant particularly for the testable PP(2v), are the repulsive particle-hole spin-isospin 
interaction and the attractive particle-particle interaction, which clearly play a decisive role in the 
concentration of the (3~ strength in the giant GT resonance, and the relative suppression of the (3 + 
strength and its concentration at low excitation energies. Together, these two ingredients are able to 
explain the values of M"q t provided the empirical parameter g pp , the strength of the particle-particle 
interaction is adjusted (but its adjusted value is, reassuringly near its expected value). Yet, the 
QRPA is often criticized for two "undesirable" features. One is the sensitivity of Mq T to the c/ pp 
value which decreases the predictive power of the method. The other is the fact that for a realistic 
value of g pp the QRPA solutions are close to their critical value (so called collapse, beyond which 
the solutions of the QRPA equations do not exist). The collapse indicates a phase transition, i.e., 
a rearrangement of the nuclear ground state. QRPA is meant to describe small deviations from the 
unperturbed ground state, and thus is not fully applicable near the point of collapse. Numerous 
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approaches have been made to extend the range of validity of QRPA, usually involving corrections 
to the quasi-boson approximation. (See, e.g., Ref. (p23|)). 

While extensive work has been done on this aspect of QRPA, the unresolved question is what 
effect on the (3(3 nuclear matrix elements, and in particular on the (3(3{Qv) ones, have the more 
complicated configurations that are not included in the QRPA (i.e., the two-quasiparticle states and 
their iterations). There is also a lack of detailed nuclear spectroscopy predictions (beyond the (3 
strength function). 

Ideally, the nuclear shell model, NSM, is the method of choice for the evaluation of (3(3 nuclear 
matrix elements. In it, one chooses a set (limited basically by the capability of the present-day 
computers) of valence single particle states. Then one finds an effective hamiltonian, based usually 
on the free nucleon-nucleon interaction, but modified to describe the effective nucleon interaction 
for that particular set. All configurations (or at least a convergent set of them) are used in the 
diagonalization of the hamiltonian, and in the evaluation of the (3(3 nuclear matrix elements. The 
method is tested, and the hamiltonian is adjusted, by requiring that it describes the spectroscopy 
(level energies and transition probabilities) of the relevant nuclei. 

Despite the tremendous advances of the computational techniques, only a limited set of single- 
particle states can be included in the NSM. The effects of single-particle states that are excluded is 
usually simulated in the NSM by using effective operators, or more simply, by using effective charges. 
Unfortunately, it is not clear what effective charges, if any (or indeed what effective operators), should 
be used in the NSM evaluation of the (3(3(§v) matrix elements Mp and Mq^, and it is not clear what 
related phenomena one can use to determine them. 

Thus, each approach, QRPA and NSM, has its strengths and weaknesses, and naturally its critics 
and defenders. Given this situation, it is customary, although not really justified, to consider the 
spread of the theoretically calculated (3(3{Qv) nuclear matrix elements as a measure of their uncer- 
tainty. Clearly, a breakthrough in the evaluation of these matrix elements, or at least in the estimate 
of their uncertainty, would be very welcome. Given the experimental effort described in the following 
sections and the importance of the problem, we hope that a comparable theoretical effort will emerge 
and result in radical improvement in the nuclear matrix element evaluation. 

To better appreciate the difficulty in evaluating the (3(3(0^) nuclear matrix elements, let us stress 
once more that the presence of the neutrino propagator leads to the appearance of the "neutrino 



potential" H(r,A), Equation |TJ. Then, 



M^ T = (f\J2^r^r+H(r lk ,A)\i) , (32) 



Ik 



M°/ = (f\J2^r+H(r lk ,A)\i) . (33) 

Ik 



Here the I, k summation is over all pairs of neutrons (or protons). Note that, due to the presence of 
H(r,A), the Fermi matrix element M F V is nonvanishing even if isospin is conserved. 

One can now expand the potential H (r, A) in multipoles corresponding to the various angular 
momenta of the intermediate odd-odd nucleus. One finds then, as expected due to the high excita- 
tion energy (or high value of the momentum q of the virtual neutrino), that many multipoles give 
comparable contributions. Moreover, the 1 + multipole, which is the only one contributing to the 
(3(3{2v), is suppressed and contributes very little. Thus, a correct reproduction of measured M"q T is 
a necessary but insufficient condition for equally successful evaluation of the M^ T and M F V nuclear 
matrix elements. 
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It has been often argued that, unlike Mq^, there is no suppression in M% T and Mjjf, and hence 
their values are less sensitive to nuclear structure details. This argument is based on the multipole 
decomposition, discussed above, in the various angular momenta and parities of the virtual states in 
the intermediate odd-odd nucleus. It turns out that the contributions of most of the mutipoles have 
the same sign, and hence do not interfere with each other. 

However, it is possible to expand the corresponding expression in an equivalent representation in 
terms of the angular momenta and parities of the pair of neutrons that are transformed into the pairs 
of protons. In this, equally valid representation, the dominant contribution of the J = + pairs is 
to a large extent cancelled by the contribution of all other, J / + , pairs which have an opposite 
sign. Thus, also in the (3/3(Qv) case the nuclear matrix elements depend on the small, and presumably 
poorly determined, pieces of the nuclear wave function. 

Despite our reservations expressed above, in Table || we compare the /3(3(0v) half-lives evaluated for 
(m u ) = 50 meV with nuclear matrix elements evaluated in the quoted references chosen to represent 
the vast literature on the subject. The spread of the calculated values for the given parent nucleus 
gives some indication of the role played by the nuclear matrix elements. On the other hand, the spread 
of half-lives along the columns in Table ^ reflects both effects of the phase space, and the nuclear 
matrix elements. The methods used to evaluate them are: truncated shell model (^); nuclear shell 
model (|80|); QRPA with the schematic 5 force interaction (a[ = 390 MeV fm 3 , recalculated for 
= 1.25) (H); QRPA with G-matrix based interaction (||); renormalized QRPA dJ,||); and QRPA 



without the p — n pairing (84). Clearly, the calculated half-life uncertainty of about an order of 
magnitude, corresponding to a factor of ~ 3 in (m u ), can be seen. 



4 P/3(0v) EXPERIMENTAL OVERVIEW AND PAST P/3(0v) EXPERIMENTS 

The various detection schemes for (3/3 are effectively outlined in Ref. ([!]) so we will only mention 
the salient points here. Over the past 15 years, the success of background reduction has resulted 
in a large number of (3(3(2v) half-life measurements. With this success behind us, the community 
is currently focusing on the more exciting goal of /3/3(Oi/). The various modes of (3(3 are separated 
by the differences in their electron sum energy spectrum and because (3(3{Qv) is identified by its 
distinguishing sum-energy peak, direct counting experiments with sufficient energy resolution are the 
focus of today's researchers. In most of these experiments, the source also serves duty as the detector. 

Geochemical and radiochemical experiments which were a mainstay of (3(3 physics through the 1970s 
and 1980s, do not discern between the different modes. Thus there is little interest in pursuing these 
techniques further. Instead relatively new technologies such as bolometers and scintillating crystals 
are receiving attention. Tracking and foil-scintillator sandwich experiments are also being pursued 
although the source and detector are separate. Amazingly, the long-time workhorse of germanium 
detectors dominates the present (3(3(Qv) decay results and comprise some of the most promising future 
proposals. 

The study of double beta decay is about suppressing backgrounds. Therefore, in this section we 
summarize the criteria that make a good (3(3 experiment and then discuss the background issues in 
general. Finally, we describe the various past and current experiments. 



4-1 Experimental Criteria 

For the best sensitivity to (m u ), one must build a detector that maximizes the (3(3{Qv) count rate while 
minimizing the background. The signal sensitivity is approximately the statistical precision (i.e., 
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square root) of the background determination. Since the number of background counts increases 
linearly with time, the decay rate sensitivity scales as the square root of time. In turn, the (m„) 
sensitivity scales as the square root of the decay rate, and therefore as the fourth root of the counting 
time. In an experimet with zero background on the other hand, the {m v ) limit scales more quickly 
as the square root of the counting time. Explicitly, the limit on (m u ) can be expressed in terms of 
experimental parameters as rt8q) 



K) = (2.50 x 10" 8 eV) 
(m„) = (2.67 x 10" 8 eV) 



W 



[fxeG^\M 0u \ 2 



1/2 



bAE 
MT 



1/4 

background limited 



W I 1 / 2 1 

x ; zero background (34) 



fxeG°»\M 0u \ 



where W is the molecular weight of the source material, / is the isotopic abundance, x is the number 
of (5(5 atoms per molecule, e is the detector efficiency, b is the number of background counts per 
kg-year-keV, AE is the energy window for (5(5{Qv) in keV, M is the mass of isotope in kg, T is 
the live time of the experiment in years, and \Mq v \ is a shorthand for the (5(5{Qv) matrix elements 



given by Equation 18. Some of the criteria that need consideration when optimizing the design of 
a (5(5{Qv) decay experiment are obvious whereas others are more subtle. It is clear from Equation 
|34| that one needs a large source mass. To reach the 50 meV region of interest indicated by the 
oscillation results, approximately a ton of isotope will be required. Also obvious is that one needs 
a reliable detector technology that preferably is easy to operate. Since the experiments are usually 
conducted in underground laboratories far from the experimenter's home, it is a great convenience if 
the experiment needs minimal maintenance. 

The search for (5(5{Qv) decay is a search for a peak superimposed on a continuum. Therefore good 
energy resolution is a must. Not only does it improve the signal to background in the peak search, 
but poor resolution would result in the (5(5(2v) tail extending up into the peak region to become a 
background itself. 

Natural radioactivity is present in all materials at some level. Thus the source and detector must 
be very low in such impurities. Furthermore, the detector must be shielded from the environment 
and its associated radioactivity. This shielding must be also radiopure. Since the total activity of 
an impure material will scale as its volume, it is usually an advantage to minimize the detector size. 
This can be most readily accomplished by employing a detector that also plays the role of the source. 

Cosmogenic activities build up in materials through nuclear reactions of cosmic ray muons and 
their secondary products, especially neutrons. These can be a significant background contribution 
both for the source and for the shielding material. Some materials have no long-lived isotopes and 
thus have a built-in safeguard against cosmogenics. For experiments that must fight this problem, 
fabricating the apparatus underground and storing materials underground can greatly reduce this 
background. Unlike solid sources, a gaseous or liquid source can be continuously purified of such 
impurities. 

Choosing an isotope with a large Q-value and matrix element improves the (m„) sensitivity for a 
given measured half-life. But the nuclear theory of some isotopes is better understood than others. 
Since the only feasible experiment sensitive to {m u ) is double beta decay and since the half-life for 



this process depends on both parameters (Equation 18), it is advantageous to use a source isotope 
for which there is confidence in the theoretical calculations. Also some nuclei, 100 Mo for example, 
have relatively fast (5(5{2v) rates with respect to the theoretically anticipated (5(5(0i/) rate for a given 
(m u ). As (5(5{2v) is a potential background, the source choice may be important for detectors with 
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modest or marginal energy resolution. 

Radiochemical and geochemical experiments operate by detecting the daughter of the decay. How- 
ever, since these techinques integrate over an exposure time, they can not identify the mode of decay. 
But if an experiment could identify the daughter in coincidence with a real-time measurement of the 
decay energy, it would have a powerful tool for rejecting many backgrounds. In this case, only (5(5{2v) 
would be a background. There are several possibilities for detecting the daughter. If the decay is to 
an excited state, one might be able to observe 7 rays that identify the daughter. (5(5 candidates are 
initially in + states and in many cases transitions to excited 2 + or + states are possible. However, 
the Q-value for these excited-state transitions is much smaller than for the ground-state transition 
and therefore the decay rate for a given {m v ) value is much less. For the 2 + case, the matrix elements 
are also much smaller due to the forbidden nature of the transition. One interesting possibility is that 
of 150 Nd where the excited state is relatively low in energy. But this excited state decays via internal 
conversion requiring the detection of a 30-keV x ray in order to observe the daughter. The most 
enticing situation however, is that of the Xe-Ba system. The optical detection of the Ba daughter ion 
might be possible. This possibility is discussed in Section 5. 



4-2 Backgrounds 

Any extraneous energy deposit in the detector near the PP(0v) Q-value will limit the sensitivity 
to (m u ). Thus any radioactive isotope with a Q-value greater than the (3(3(0u) endpoint may be a 
potential background. Since the number of radioactive isotopes decreases with increasing Q, it is 
desirable to select a (5(5 candidate with as large a Q-value as possible. (5- and a-emitting decays 
are easy to shield and thus are only a problem if they occur within the detector or on its surface. 
Penetrating 7 rays pose a more difficult problem. In this subsection, we consider the origin of various 
sources of background and some techniques for their mitigation. 



4.2.1 Natural Activity 

The naturally occurring isotopes of U and Th and their daughters are present as impurities in all 
materials at some level. The half-life of the chain patriarch is comparable to the age of the universe 
but very short compared to the half-life sensitivity of the /3/3(0^) experiments. Therefore even a small 
quantity of U or Th will create a significant background. In particular, 214 Bi and 208 T1 have large 
Q-values and the decay spectra will overlap the endpoint of almost all the f5(5(Oi / ) candidates. Even 
tracking experiments have difficulty with these two isotopes due to their (5 decays which are promptly 
followed by internal conversion resulting in a two-electron event that mimics (5(5(Qv) or (5(5{2v). 
Careful selection of materials and purification have been successful, if difficult, solutions to these 
problems. In the past decade, great strides have been made in purifying some materials. Of particular 
note are liquid scintillator (86), electroformed Cu(87), and CVD nickel (p^). The radioactive chains 



may or may not be in equilibrium depending on the sample's history because chemical treatment or 
purification can disproportionately eliminate the daughters. 

Radon gas, either 222 Rn or 220 Rn, is especially intrusive and may infiltrate into a detector's sensitive 
region. These parents to 214 Bi and 208 T1 are mobile and diffuse through many materials. Their 
daughters tend to be charged and stick to dust or any other electrostatic surfaces. Many experiments 
eliminate Rn from the detector vicinity by purging the volume immediately surrounding it with N2 gas 
that has boiled-off from a liquid nitrogen (LN) supply. Because Rn freezes out at LN temperatures, 
the boil-off gas tends to be very low in Rn especially compared to the laboratory air being displaced. 
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Some groups have also installed charcoal Rn scrubbers into the laboratory airstream. 

Activities such as 3 H, 14 C, and 40 K are also naturally present but their Q-values are too low to 
interfere with f3(3(0v) experiments. However, many Pf3(0u) experiments double as dark matter ex- 
periments by studying the low end of their energy spectrum for possible elastic scattering of WIMPS. 
These experiments have to consider a larger pool of potential background isotopes. 

4.2.2 Cosmogenic and Induced Activities 

Long-lived radioactive isotopes can be produced through various nuclear reactions. Many of these 
isotopes can have decay energies that exceed the (3(3(0v) Q-value and thus can create a background. 
Such activities can be produced in the detector or the shielding material. In particular, long-lived 
isotopes of the source element can be very difficult to remove by purification. The troublesome 
68 Ge 271-day activity that builds up in Ge detectors is a good example of this problem. Short-lived 
activities can also create background if created in situ while the experiment is operating. There 
are several nuclear processes that should be considered and we discuss them in this section. The 
magnitude of the background due to each process is material dependent, and the flux of the projectiles 
that induce the activities depends on the environment. In fact, many processes are greatly minimized 
by going underground where the cosmic ray flux is decreased. 

Neutron capture produces 7 rays and frequently a radioactive isotope. Since neutrons are neutral 
and difficult to identify with anti-coincidence detectors, they can be a significant problem. On the 
Earth's surface most neutrons arise from the cosmic-ray hadronic component. In shallow laboratories, 
secondary neutrons from cosmic-ray muon interactions can form a large contribution to the total 
neutron flux. In deep sites, however, the neutron flux is dominated by (a, n) reactions and fission 
neutrons from the laboratory's rock walls. Siting deep underground, covering the walls with shielding 
material to reduce the overall flux inside the laboratory, and placing neutron shielding around the 
detector can help control this background. 

Fast neutron reactions also need consideration. 68 Ge for example is produced by fast neutron 
( 25 MeV) interactions on stable Ge isotopes. Above ground the dominant source of these fast 
neutrons are secondaries produced by cosmic rays. Once the material is taken deep underground, 
the problem is mitigated for the most part, as only the residual surface production remains. But to 
estimate the underground production rate at a given depths, one requires neutron flux data. There are 
several integral measurements of the neutron flux underground along with estimates of the neutron 
production due to muons. However, since the neutron flux falls very quickly with energy, it is difficult 
to deduce the higher energy flux from these measurements. A summary of the past measurements 
and a calculation of the fast neutron flux is given in Ref. (|89| ) and a program to measure the flux is 
being developed (|90|) . 

Muons and muon induced electromagnetic showers can produce background also. Going deep 
underground reduces the flux, and veto systems surrounding the detector are used to eliminate any 
prompt activity observed in coincidence. However, inelastic fi scattering and \x~ capture produce 
delayed events inside the detector after the \x signal. If this time delay is too long or if the fi flux is 
too high, anti-coincidence techniques won't suffice. Additional depth can remedy this problem. For 
high-Z materials [i~ capture dominates over /i decay ( |9l|) and the neutron multiplicity is of order 
1, with an energy spectrum extending to many tens of MeV (|92"|). Many of the solar neutrino, dark 
matter, and double beta decay experiments have analyzed the possible spallation products that might 
be produced. (See for example Ref. ( |93| ) which gives a brief general discussion of the topic.) 
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4.2.3 Artificially produced activity 

Artificial radioactive isotopes can also be present in materials. For example ~ 10 15 Bq of 239 > 240 p u 
from the above-ground testing of nuclear weapons coats the surface of the Earth. Nuclear accidents 
such as that at Chernobyl have introduced long-lived isotopes e.g., 137 Cs, 90 Sr, and Pu into the 
environment also. When considering what backgrounds might be present, it is thus prudent to 
consider these exotic possibilities. The noble gas radioactive isotopes 42 Ar and 85 Kr arise from the 
venting of reactors and atmospheric testing; these must be considered by Xe f3/3 experiments. 



4.2.4 (3f3(2v) as a background 

When searching for the /3/3(Oi/) peak, one must consider (3(3(2v) decay as a potential ultimate back- 
ground. Near the endpoint energy (Q), the /3/3(2v) spectrum has very little strength. But since Ty 2 
is much shorter than T®%-> the effect of resolution (shown in Fig. |]) must be considered. Roughly 
speaking, the f3/3(2v) counts within one peak width (AE) centered on Q will contribute to the (3f3(Qv) 
peak region and be a background. The fraction (F) of the [3[3{2u) counts in the peak region can be 
approximated by 

F=I«1 (35) 
m e 

where 5 (= AE/Q) is the FWHM energy resolution expressed as a fraction and m e is the electron 
mass. The coefficient 7 is for a resolution of 5%. This coefficient depends moderately on resolution 
and is 8.5 (5) at 1% (10%). An expression for the /3/3(0i/) signal (S) to (3(3{2v) background (B) ratio 
can then be written 

S_ m e r 0l , _ m e Ty2 

B~ 7Q5*T 2u ~ 7QS e T°J 2 [ j 

Although this approximation cannot replace a Monte Carlo simulation of an experiment's perfor- 
mance, it clearly indicates that good energy resolution is critical. But, in addition, the ratio of Ty 2 
(for (m„) = 1 eV) to Ty 2 can vary from 5,000 to 100,000 depending on the isotope's Q-value and 
matrix element estimate. Therefore, the choice of isotope needs consideration. 

For a S/B of 1, the (m„) sensitivity limit due to the /3(3(2v) background can be estimated as 

, , 2 7Q6« G 2u \M 2u \* 

{mv) ~ —gZ W^- (37) 

Note that using an asymmetric window defined by Q < E < Q + AE/2 reduces S by a factor of 
2 but decreases B by a factor of ~16. This is exploited by some experiments where the resolution is 
not ideal. A previous calculation of the f3[3{2v) contribution to the upper half of the /3f3(0u) window 
differs from this estimate. The author of Ref. ( |94| ) agrees with the present result. 



4-3 The Past Experiments 

We have written this section to embellish the work of Ref. (|T|), not replace it. Thus we have 
de-emphasized material already covered there. There has been impressive progress in systematically 
cataloging j3j3{2v) rates as detailed in Ref. (|28| ) and summarized in Table |l|. In this summary table, we 
have listed average values of the 10 measured /3/3(2i/) half-lives and deduced \M 2v \ values from them. 
The quoted value for each parent nucleus is the weighted average of the chosen measurements. We 
included in the average selected measurements with quoted uncertainties small enough to significantly 
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affect the average. To assign the individual uncertainty associated with each measurement, we first 
separately averaged the asymmetric statistical and systematic errors, and then added the two classes of 
errors in quadrature. The summarized Ty 2 are quoted at the 68% confidence level. In the individual 
measurements where the uncertainty range was quoted at a different confidence level, we scaled 
the uncertainties so they would correspond to a 68% confidence range. In the case of 96 Zr the 
measurements are inconsistent and we chose the spread of the measurements as an indication of the 
uncertainty in the measured Ty 2 . In the case of 100 Mo, one measurement was very different than 
the others and we did not use it in the average. The nuclear matrix elements were deduced using the 
phase-space factors of Ref . ( p2|) . This procedure is somewhat arbritrary but the details on the half-life 



measurements can be found in the quoted references. Although 136 Xe has not had its Ty 2 measured 
yet, it is an isotope under consideration for some of the future big experiments and therefore its Ty 2 
is of importance. We list its limit in the table. 

It is worthwhile to note that in one case, for the (5(5 decay of 100 Mo, the transition to the excited 
+ state at 1.13 MeV in 100 Ru has been observed (95,96|). That state deexites by emission of two 7 
rays; their observation serves as a convenient and clean signature of the decay. The resulting averaged 
half-life, l\/ 2 = (6.8 ± 1.2) x 10 20 y, corresponds to a matrix element of similar magnitude as the 2v 
transition to the 100 Ru ground state, in agreement with expectations (|97|). Similar transitions are 
possible in other (5(5 decay candidates, but have not yet been observed. 

The (5(5(0v) half-life limits have also improved for many isotopes and they are summarized in Table 
||. The spread of calculated T®Y 2 given in Table § gives an indication of the uncertainty in (m u ) due 
to the uncertainty in the nuclear physics. As mentioned earlier the spread in the half-lives is about 
an order of magnitude and thus the spread in deduced (m u ) would be about a factor of 3. The best 
limits come from the 76 Ge experiments and they indicate that {m u ) < 0.3 - 1 eV. 

The half-life limits for the Majoron mode are also improving. Like the (5(5(2v) mode, the /5(5(0is, x) 
spectrum is also a continuum and therefore the limits, summarized in Table ||, are more comparable to 
(5[5{2v) than to /3/3(Oi^) . We have only considered Majoron decay modes that emit a single Majoron, 
and list in Table || the half-life limits and the corresponding (^ x )values. (See Equation 31 for the 
relationship between T®J 2 X and {g v ,x)-) Note that, naturally, the deduced (g u ,x) limits depend on the 
nuclear matrix elements. 

Interestingly, the best constraint on (g u ,x) comes from 128 Te which has the longest measured total 
half-life. Furthermore, its Q-value is very low, resulting in a relative enhancement of the phase space 
factors for the (5(5(01;, x) and f5(5(0v) modes compared to the (5(5(2v) mode. Therefore, even though 
the observed rate is most likely due to the /3/?(2i/)decay, a conservative assumption is to assign all the 
rate to an exotic mode when estimating parameter limits. In the case of 128 Te, the limit on {g U)X ) of 
< 10 -5 is the best, and for (5(5(§v) the limit of ~ 1 eV is competitive. 

In the remainder of this section we discuss selected recent (5(5 experiments that can be considered 
effective prototypes for future programs. We have ordered the discussion in the current section to 
parallel Section 5 according to experimental technique. 



4.3.1 MIBETA 

The MIBETA experiment ( |101[ ) used Te02 crystals as bolometers. These detectors exploit the low 
heat capacity of the crystals at low temperature. A small energy deposit therefore results in a 
significant temperature increase of the crystal. The experiment consisted of an array of 20 crystals 
totaling 6.8 kg. Since 130 Te is 33.8% naturally abundant, enrichment was not necessary although 
two crystals were enriched to 93% in 130 Te and two others were enriched to 95% in 128 Te. The 
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crystals were arranged into a tower of 5 layers of 4 detectors within a dilution refrigerator 3500 mwe 
underground at the Laboratori Nazionali del Gran Sasso. The tower frame was made of OFHC copper 
with the crystal supports composed of Teflon. The temperature sensors were neutron transmutation 
doped germanium thermistors. Old Roman lead (<4 mBq/kg 210 Pb) was placed inside the cryostat 
surrounding the tower. The dilution refrigerator itself was shielded with low-activity lead (16 ± 4 
Bq/kg 210 Pb). The array was operated at a temperature of ~12 mK with an array-averaged resolution 
of «8 keV FWHM (0.3%) at the /3f3(0v) endpoint energy of 2.529 MeV. Since thermal detectors are 
sensitive over their entire volume, they are susceptible to surface contamination and indeed these 
crystals did observe a surface a activity that contributed significantly to the (3(3{Qv) window. The 
cosmogenically produced activities in Te are short-lived and therefore posed no significant problem. 



4.3.2 Gotthard tunnel 

The Gotthard Xe experiment ( |102[ ) used a 5-atm gas time projection chamber with 3.3 kg of 62.5% 
enriched 136 Xe. The tracking feature of the detector permitted the identification of two-electron 
tracks indicative of f3f3 decay. The energy resolution at the /3/3(0^) endpoint (2.481 MeV) was « 165 
keV FWHM (6.6%). The dominant background for /3/3(Oz/) was concluded to be Compton scattered 
electrons from natural 7 activities. These electrons were occasionally misidentified as two-electron 
events. Cosmogenic activities are not a serious issue for Xe experiments, because there are no long 
lived Xe isotopes and liquid or gaseous Xe can be continuously purified of non-Xe isotopes. 



4.3.3 Heidelberg-Moscow and IGEX 



The Heidelberg-Moscow ( 104 , 57|) (hereafter referred to as HM) and IGEX(|99|) (International Germa- 



nium Experiment) collaborations both used Ge detectors 86% enriched in 76 Ge. HM used 125.5 moles 
of active material whereas IGEX used 90 moles. With comparable masses and run times, the results 
from the 2 experiments are similar with HM posting a modestly better limit. Since PP(0u) 
events produce localized ionization in the detectors and many backgrounds (e.g., Compton scattering 
of 7 rays) produce multi-site energy deposits, both experiments used pulse-shape discrimination to 
reduce background. 

HM identifies the radioactivities that contribute to the data by their associated peaks within the 
spectrum. With these identifications, the response of the detectors to was simulated with Monte 
Carlo. A fit to the actual data using the Monte Carlo spectra provided an indication as to the 
location of the activity. The conclusion of that work, described in Ref. (57), is that the copper parts 
of the cryostat contained the majority of the background sources. 

IGEX performed a measurement of the cosmogenic activity produced in Ge crystals and compared 
that to a calculation of the expected rate based on the measured surface neutron flux and neutron 



interaction cross sections (105). The calculated rates and measured rates agreed well. It was found 
that initially 68 Ge was the dominant cosmogenic activity with the longer-lived 60 Co dominating at 
later times. The rate of 68 Ge can be determined within each crystal by measuring the intensity of 
the 10.4-keV x-ray peak. It was found that the count rate in the (3(3{Qv) window could be mostly 
attributed to radon intrusion. However, the authors concluded that with a reduction of radon, Ge 
activation isotopes would be the limiting background source in the IGEX experiment. 

The two experiments quote similar background levels in the (3(3{Qv) region of ~0.20 counts/(keV-kg-yr) 
before pulse shape discrimination and ~0.06 counts/ (keV-kg-yr) after. However, remarkably, the two 
collaborations have come to very different conclusions as to the composition of the limiting component 
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of the background in these experiments. Resolving this issue is one of the most critical debates in 
experimental 0/3(01/) research today as the design of the next generation of Ge detector experiments 
depends heavily on its outcome. (See Section 5.1.) 

During the final preparations of this manuscript (Jan. 2002), the paper Ref. ( |106| ) appreared in 
print. It uses the HM data to claim evidence for 00(Ou) in 76 Ge with a iffa = (0.8 - 18.3) x 10 25 
y. If true, this result would be extremely important and hence requires extensive substantiation and 
review. However, by itself, the paper does not sufficiently support the claim (|107|) . Such deficiences 
are not necessarily indications that the claim is wrong, but they indicate that the assessment of this 
result by the 00 community will be some time in coming. In particular, the questions raised in ( |107| ) 
should be answered first by the authors of Ref. (|106| ) . 



4.3.4 UCI, ELEGANTS and NEMO 

The University of California at Irvine time projection experiment (TPC) measured several isotopes. 
Each sample was a few tens of grams and placed as a thin foil upon the central electrode of a TPC. 
On either side of this source plane were drift regions for recording a three-dimensional image of the 
ionization trails produced by the 00(2u) electrons. The experiment was very successful utilizing the 
tracking capability to determine several kinematic parameters characterizing events. This information 



was critical to reducing backgrounds to a level that allowed the first direct detection of 00 (|6C, 108 , 
|109| ) . The drawback of the design is the limited amount of source mass as compared to the size and 
complexity of the detector. 

The best limits on 100 Mo 00(Ov) decay come from the ELEGANTS (ELEctron Gamma-ray Neu- 
trino Telescope) experiments (|100[ ). The emitted electrons in this experiment traversed drift chambers 
for measuring their trajectories and then passed into plastic scintillator in order to measure their en- 
ergies and arrival times. Nal arrays surrounded the apparatus to provide 7- and x-ray observation. 
Copper and lead shielding enclosed the detectors. The 171-g 100 Mo source was two thin foils situated 
between the drift chambers. This detector had a diminished 00(Ov) detection efficiency compared to 
the Ge detectors or bolometers. However, it had additional background rejection power because of 
its measurement of several kinematic parameters. The dominant backgrounds in the 00(§v) window 
were identified to be 214 Bi and 208 T1 contained in the source film and detector elements. 



The NEMO-2 experiment (Neutrino Ettore Majorana Observatory) (|61,62, 72, 110 , |65| ) has analyzed 



the 00(2v) rate for several isotopes. The detector had a tracking volume of 1 m 3 of He gas with two 
sides covered by scintillator calorimeters. The tracking volume was bisected by a thin source plane 
and consisted of frames containing crossed Geiger cells. An electron was defined by a track passing 
from the source foil to the calorimeter. The three dimensional track measurements were made using 
the drift times and plasma propagation times of the Geiger cells. The energy was determined with 
the calorimeters. The various source foils weighed up to about 175 g. 

These tracking experiments all had small source masses, modest energy resolution, and a complex 
apparatus. As a result, it will be a challenge for future efforts modeled on these designs to be 
competitive in the search for /?/3(0^)decay. 



4.3.5 Scintillating crystals 

There has been some progress making large scintillating crystals with an appreciable amount of 00 
isotope contained. The Bejing group used CaF2 (|^) to study 48 Ca and placed a lower limit on the 
00(Ou) decay rate. A Kiev-Firenze collaboration(|73D has used 116 CdW04 scintillators to measure the 
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(3(3(2u) half-life of 116 Cd and placed limits on the (3(3(Qv) and /3f3(0v, x) modes. These experiments are 
important in that future experiments are being planned that will exploit similar crystal technologies. 



5 FUTURE P/3(0v) EXPERIMENTS AND PROPOSALS 

5.1 The Various Proposals 

It has not been possible to realize an experimental concept that features all the criteria described in 
Section 4. These criteria are frequently incompatible and thus no past experiment or future proposal 
has been able to optimize each simultaneously. We are aware of over 10 ideas or proposals for 
(3j3{Qv) experiments summarized in Table ||. The fourteen listed proposals are divided into 2 groups 
according to source mass and arranged alphabetically within each group. Each proposal chooses a 
different approach to attain this optimization. Five of them are substantially developed and have 
the potential to reach the crucial 50 meV region. Although we will briefly mention all proposals, we 
devote particular attention to these five: CUORE, EXO, GENIUS, Majorana, and MOON. 



5.1.1 CUORE 

The success of the MIBETA experiment has resulted in the CUORE (Cryogenic Underground Ob- 



servatory for Rare Events) proposal ( [1251) . 1000 Te02 crystals of 750 g each would be operated as a 
collection of bolometers. The detectors will be collected into 25 separate towers of 40 crystals. Each 
tower will have 10 planes of 4 crystals each. One such plane has already been successfully tested and 
a single tower prototype referred to as CUORICINO has been approved. 

The energy resolution at the P(3(0v) peak (2.529 MeV) is expected to be about 5 keV FWHM (« 
0.2%). A low energy threshold of 5-10 keV is anticipated, and thus the experiment will also search 
for dark matter. The background has been measured in the first plane to be ~0.5 counts/(keV-kg-y). 
However a major component of this background is due to a surface contamination arising from the use 
of cerium oxide polishing compound which tends to be high in thorium. With this problem solved, 
the experimenters project a conservative estimate of the background to be ~0.01 counts/(keV-kg-y). 

A major advantage of this proposal is that the natural abundance of 130 Te is 34% and, thus, 
no enrichment is needed resulting in significant cost savings. As with MIBETA, the cosmogenic 
activities within the Te02 crystals are not a serious concern. On the other hand, the crystal mounts 
and cryostat form a significant amount of material close to the bolometers. Much of the cryostat is 
shielded with Roman period lead but a fair quantity of copper and Teflon remain close to the crystals. 



5.1.2 EXO 

The Enriched Xenon Observatory (EXO) (|126|) proposes to use up to 10 t of 60-80% enriched 136 Xe. 
The unique aspect of this proposal is the plan to detect the 136 Ba daughter ion correlated with the 
decay. If the technique is perfected, it would eliminate all background except that associated with 
/3(3(2v). The real-time optical detection of the daughter Ba ion, initially suggested in @), might 
be possible if the ion can be localized and probed with lasers. The spectroscopy has been used for 
Ba + ions in atom traps. However, the additional technology to detect single Ba ions in a condensed 
medium or to extract single Ba ions from a condensed medium and trap them must be demonstrated 
for this application. To optically detect the alkali-like Ba + ion, it is excited from a 6 2 S 1 / 2 ground 
state to a 6 2 Pi/2 with a 493-nm laser. Since this excited state has a 30% branching ratio to a 5 4 D 3 / 2 
metastable state, the ion is detected by re-exciting this metastable state to the 6P state via a 650-nm 
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laser and then observing the resulting decay back to the ground state. This procedure can be repeated 
millions of times per second on a single ion and produce a significant signal. 

EXO is presently considering two detector concepts: high-pressure-gas Xe TPC or Liquid Xe (LXe) 
scintillator. The TPC baseline design consists of two 35 m 3 modules at ~20 atm for a total of 8.4 t of 
Xe. The Xe would be contained in a non-structural bag within a pressurized buffer gas to constrain 
the Xe to the active region. The spatial resolution and typical (3 particle range (5 cm) will permit the 
identification of the high-ionization-density points at the terminus of the beta tracks, aiding in the 
separation of two-electron events from one-electron backgrounds such as Compton scatters. Upon a 
trigger of an event near the /3/3(Qv) peak energy, 2.481 MeV, the lasers are directed to the decay point 
to excite the Ba + ion. One complication is that the P(3(0v) of 136 Xe produces a Ba ++ ion whereas 
the spectroscopy requires a Ba + ion. Because Xe is a tightly bound atom, charge exchange with the 
Ba ion is unlikely and a quenching gas is required to neutralize one stage of ionization. 

The EXO LXe concept has the advantage of being much smaller than the TPC due to the high 
density of LXe. The scintillation readout has better energy resolution but cannot spatially resolve the 
high ionization points. The higher density makes the scattering of the laser light too great to optically 
detect the Ba + in-situ. However, once the Ba ion is localized via its scintillation and ionization, it 
might be extracted via a cold finger electrode coated in frozen Xe (M. Vient, unpublished observation, 
1991). The ion is electrostatically attracted to the cold finger which later can be heated to evaporate 
the Xe and release the Ba ion into a radio frequency quadrupole trap. At that point, the Ba ++ is 
neutralized to Ba + , laser cooled and optically detected. The efficiency of the tagging has yet to be 
demonstrated and is a focus of current research. 

The collaboration is currently performing experiments to optimize the energy resolution for both 
configurations. The resolution is a critical parameter as f3[3{2v) would then be the lone background 
if the Ba tagging is successful. Tests to determine the viability of the Ba extraction process are also 
being performed. The EXO collaboration has received funding to proceed with a 100-kg enriched Xe 
detector without Ba tagging. This initial prototype will operate at the Waste Isolation Pilot Plant 
(WIPP) in southern New Mexico. 

5.1.3 GENIUS 

The progress and understanding of Ge detectors has been developed over more than 30 years of 
experience. The potential of these detectors lie in their great energy resolution, ease of operation, and 
the extensive body of experience relating to the reduction of backgrounds. This potential is not yet 
exhausted as is evidenced by the GENIUS and Major ana proposals that build on the experimenters' 
previous efforts. 

The GENIUS (GErmanium Nitrogen Underground Setup) (|11S| ) proposal has evolved from the 
Heidelberg-Moscow (HM) experiment. The driving design principle behind this proposed Ge detector 
array experiment is the evidence that the dominant background in the HM experiment was due to 
radioactivity external to the Ge. (The reader should contrast this with the motivation for the design 
of the Majorana proposal described below.) An array of 2.5-kg, p-type Ge crystals would be operated 
"naked" within a large liquid nitrogen (LN) bath. By using naked crystals, the external activity 
would be moved to outside the LN region. P-type crystals have a dead layer on the external surface 
that reduces their sensitivity to external and a activity. Due to its low stopping power, roughly 12 
m of LN is required to shield the crystals from the ambient 7-ray flux at the intended experimental 
site at Gran Sasso. By immersion in LN, the optimal operating temperature is maintained without 
a bulky cryostat and a test of the naked operation of a crystal in a 50 1 dewar has been successful 
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( |127| ), The results indicate that the performance of the detector was comparable to those operated in 
a conventional vacuum-tight cryostat system. Their measurements also indicate very little cross talk 
between naked detectors and that an extended distance (~6 m) between the FET and the crystal 
does not degrade the signal. 

The proposal anticipates an energy resolution of ~ 6 keV FWHM (0.3%) and a threshold of 11 
keV. The value of this low threshold is set by x rays from cosmogenic activities. Using 1 t of 86% 
enriched Ge detectors, the target mass is large enough for dark matter studies. In fact a 40-kg nat Ge 
proof-of-principle experiment has been approved for dark matter studies. 



5.1.4 Major ana 

The Majorana proposal ( |121|) (named in honor of Etorre Majorana) involves many of the IGEX 
collaborators. Their analysis indicated that 68 Ge contained within the Ge detectors was the limiting 
background for their /3/3(Oi/) search. (Contrast this with the GENIUS approach described above.) 
The proposal's design therefore emphasizes segmentation and pulse shape discrimination to reject 
this background. The electron capture of 68 Ge is not a significant problem but 68 Ge decays to 
the (3 + emitting 68 Ga. This isotope can create background in the (3{3(0v) window if one of the 
annihilation 7 rays converts within the crystal. The energy deposits of the positron and 7 ray may 
polute the peak window in energy, but the deposits will be separated in space. In contrast, a f3f3{Qv) 
event will have a localized energy deposit. Segmentation of the crystals permits a veto of such 
events. Furthermore, distinct ionization events will have a different pulse shape than a localized 
event. Therefore pulse shape analysis can also help reject background. Majorana plans to use 210, 
86% enriched, segmented Ge crystals for a total of 500 kg of detector. The cryostat would be formed 
from very pure electroformed Cu (< 25 //Bq/kg 226 Ra, 9 /iBq/kg 228 Th)(]87|). 



5.1.5 MOON 

The MOON (Mo Observatory Of Neutrinos) proposal (|122j ) plans to use 100 Mo as a j3j3(0u) source and 
as a target for solar neutrinos. This dual purpose and a sensitivity to low-energy supernova electron 
neutrinos (128) make it an enticing idea. 100 Mo has a high Q-value (3.034 MeV), which results in a 



large phase space factor and places the (3f3(0v) peak region well above most radioactive backgrounds. 
It also has hints of a favorable \M$ V \ but unfortunately it has a fast Ty 2 . The experiment will make 
energy and angular correlation studies of (3j3 to select /3f3(0v) events and to reject backgrounds. The 
planned MOON configuration is a supermodule of scintillator and Mo ensembles. One option is a 
module of plastic fiber scintillators with thin (0.03 g/cm 2 ) layers of claded Mo, which are arranged 
to achieve a position resolution comparable to the fiber diameter (2-3 mm). A total of 34 tons of 
natural Mo would be required. 

As a solar neutrino detector, 100 Mo has a low threshold: 168 keV, and the estimated observed 
event rate is «160/(ton 100 Moccfoiyear) without neutrino oscillations. It is sobering to realize that 
the primary background for the delayed-coincidence solar neutrino signal is accidental coincidences 
between 0f3(2v) decays. 

The project needs Mo and scintillator radioactive impurity levels of better than 1 mBq/ton. This 
can be achieved by carbonyl chemistry for Mo and plastics can be produced cleanly. However, the 
total surface area of the Mo-scintillator modules is ~26000 m 2 . Dust, being electrostatically charged, 
tends to garner Rn daughters and becomes radioactive. Keeping these surfaces clean of dust during 
production and assembly will be a challenge. Liquid scintillator and bolometer options that would 



24 



Elliott and Vogel 



avoid this large surface area are also being considered. The simulations of the design indicate that the 
energy resolution for the (3(3{Qv) peak will be ~7% which is at the upper end of the range of feasibiltiy 
for a sub 50 meV (m u ) experiment. The bolometer option would also remove the resolution concerns. 
Use of enriched 100 Mo is feasible, and would reduce the total volume of the detector and source 
ensemble resulting in a lower internal radioactivity contribution to the background by an order of 
magnitude. 



5.1.6 OTHER PROPOSALS 

There are too many proposals for detailed description so we have summarized those of which we 
are aware in Table || and mention them here. The CAMEO proposal (114) would use 1000 kg 



of scintillating 116 CdW04 crystals situated within the Borexino apparatus. The Borexino liquid 
scintillator would provide shielding from external radioactivity and light piping of crystal events to 
the photomultiplier tube (PMT) array surrounding the system. Early phases of the program would 
use the Borexino counting test facility. Similarly, the CANDLES proposal( |115D (CAlcium floride for 
study of Neutrino and Dark matter by Low Energy Spectrometer) plans to imerse CaF2 in liquid 
scintillator. The scintillation light from the (3(3 of 48 Ca will be detected via PMTs. The low isotopic 
abundance (0.187%) of 48 Ca requires a very large operating mass. Two groups (119, 12~0| ) have been 
studying the use of GSO crystals (Gd2SiOs:Ce) for the study of 160 Gd. 



COBRA (CdTe O neutrino double Beta Research Apparatus) ( |lll|) would use CdTe or CdZnTe 
semiconductors to search for (3(3(Qv) in either Cd or Te. 1600 1-cm 3 crystals would provide 10 kg of 
material. GEM is a proposal ( |H7|) that is very similar to that of GENIUS. However, much of the LN 
shielding would be replaced with high purity water. 

The Drift Chamber Beta-ray Analyzer (DCBA) proposal ( [11 2D is for a three-dimensional tracking 
chamber in a uniform magnetic field. A drift chamber inside a solenoid and cosmic-ray veto counters 
comprises the detector. Thin plates of Nd would form the source. The series of NEMO experiments 
is progressing with NEMO-3 (|113| ) beginning operation in 2002. In concept the detector is similar to 
NEMO-2. That is, it contains a source foil enclosed between tracking chambers that is itself enclosed 
within a scintillator array. NEMO-3 can contain a total of 10 kg of source and plans to operate with 
several different isotopes, but with 100 Mo being the most massive at 7 kg. The collaboration is also 
discussing the possibility of building a 100-kg experiment that would be called NEMO-4. 

There are two additional groups proposing to use 136 Xe to study (3(3(Qv). Caccianiga and Gi- 



ammarchi (123) propose to dissolve 1.56 t of enriched Xe in liquid scintillator. The XMASS (124) 
collaboration proposes to use 10 t of liquid xenon for solar neutrino studies. The detector would have 
sensitivity to (3(3(0v). 



5.1.7 ISOTOPE ENRICHMENT 

Only the enrichment facilities of Russia can enrich materials in the 1-ton quantities that are required 
for the future proposals. It has several centrifuge enrichment facilities. One of these, the Electro 
Chemical Plant (ECP) in Krasnoyarsk, can produce ~30 kg/year of enriched Ge material in appa- 
ratus that have not been used for uranium enrichment. With some modest improvements, they could 
increase production to ~200 kg/year. The 136 Xe production rate at these facilities is estimated to 
be 2 t/year. With several plants throughout Russia and fairly easy expansion, the total capacity is 
large enough that two samples of isotope could be produced simultaneously. 
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We have reviewed the motivation, present status, and future plans of the search for /3/3(Oi/) decay. 
Seeing (3/3(0v) decay would be a remarkable physics result, with important consequences for neutrino 
physics in particular, and for the hunt for 'physics beyond the Standard Model' in general. 

For N = 3 neutrino flavors and mass states, the neutrino mass matrix contains, in general, nine 
parameters. Three mixing angles, two Am 2 mass differences, and one CP violation phase can be, in 
principle, determined in neutrino oscillation experiments. Remarkable progress has been made lately 
in determining some of these parameters, and great effort is being devoted to verify the discoveries 
already made and refine and extend the search for neutrino oscillations further. 

The remaining parameters of the neutrino mass matrix, the absolute mass scale and the two 
Majorana phases, can be determined or severely constrained only by the observation of (3(3{Qv) decay, 
and/or by further progress in the tritium endpoint neutrino mass experiments. Thus, the search for 
the (3(3(0v) decay has become one of the critical issues of particle physics today. 

We should note here again that if a non-zero m u is observed in either (3(3{Qv) decay or in tritium 
beta decay (or, ideally, in both), the nuclear matrix element \Mq v \ issue will become critical. To 
dramatize this problem, consider the following possibility. Suppose a non-zero {m v ) of about 100 
meV is indicated by an upcoming 76 Ge experiment. The anticipated experimental uncertainty is 
approximately ±25 meV, and the matrix element uncertainty of a factor of 2-3 would dominate the 
total uncertainty on the mass parameter. Hence the possible range of effective neutrino mass would 
be 33 meV < (m u ) < 300 meV. Keep this range in mind and look at the right hand side of Fig. ||. One 
sees that this matrix element uncertainty alone is large enough to compromise physics conclusions 
regarding the mass heirarchy at small minimum neutrino masses. Therefore, we use this opportunity 
to express our hope that the improved sensitivity of the upcoming (3(3(0v) experiments will promote 
an increased interest in the nuclear theory of double beta decay. 

Altogether, the field of (3(3 decay has seen a great rebirth of interest stimulated by the recent 
neutrino physics results. The oscillation experiments indicate that neutrinos do have mass and, in 
particular, at least one neutrino has a mass greater than m sca i e f«50 meV. The upcoming (3(3(0u) 
experiments will have a sensitivity to (m u ) values less than this critical mass scale. This is a very 
exciting time for (3(3{Qv) research as it is reasonable to hope for a positive result within the coming 
decade. 
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Figure 1: Illustration of the spectra of the sum of the electron kinetic energies K e (Q is the endpoint) 
for the (3(3{2u) normalized to 1 (dotted curve) and (3(5{^)v) decays (solid curve). The (3(3{0i>) spectrum 
is normalized to 10~ 2 (1CP 6 in the figure inset). All spectra are convolved with an energy resolution 
of 5%, representative of several experiments. However, some experiments, notably Ge, have a much 
better energy resolution. 
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Figure 2: "Moore's law" of (3(3{Qv) decay: the limit of the effective neutrino mass vs. time. The 
corresponding experiments are denoted by the symbol for the initial nucleus. The uncertainty in the 
nuclear matrix elements is not included in this illustration. The gray band near the bottom indicates 
the neutrino mass scale J Am^ tm . 
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Figure 3: Effective mass {m v ) as a function of the smallest neutrino mass m m i n . The left panel is 
for the normal mass hierarchy, as indicated in the insert (not to scale), and the right panel is for the 
inverted hierarchy. Both panels are evaluated for the LMA solar solution with Am 



atrii 



2.4 x II) :i 

eV 2 , Am 2 ^ = 4.5 x 10~ 5 eV 2 , and \U e 2\ 2 = 0.3. The full lines show (m u ) max and {m u ) m i n , defined in 
Equation 21, for U e % = and the dashed lines use the maximum value | f/ e 3 1 2 = 0.025 allowed by the 
Chooz and Palo Verde reactor experiments fll2], p|). 
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Table 1: Summary of experimentally measured j3j3{2v) half-lives and matrix elements ( 136 Xe is an 
important exception where a limit is quoted). See text for a discussion of the averaging procedure. 
In the Isotope column, the symbol * is used for inconsistent results in which case the uncertainty 
reflects the spread in the measured values. In the references column, the ' is used to indicate the 
outliers which were not used in the averaging. The nuclear matrix elements were deduced using the 



phase-space factors of Ref. ( P2| ) using the mean Tfj 2 ■ 



Isotope 


t?; 2 (y) 


48 Ca 


(4.2 ± 1.2) x 


10 19 


76 Ge 


(1.3 ±0.1) x 


10 21 


82 Se 


(9.2 ± 1.0) x 


10 19 


96 Zrt 


x 


10 19 


100 Mo 


(8.0 ±0.6) x 


10 18 


116 Cd 


(3.2 ±0.3) x 


10 19 


128 Te (l) 


(7.2 ±0.3) x 


10 24 


130 Te (2) 


(2.7 ±0.1) x 


10 21 


136 Xe 


> 8.1 x 10 20 (90% 


CL) 


150 Ndt 


7.0^ 3 8 x 


10 18 


238tj(3) 


(2.0 ±0.6) x 


10 21 



References 



M'q T (MeV^ 1 ) 



00' 



©0) 



(0,00 



0,000,0,[7(]),(0) 



mm 
0) 



dZ3) 



(|7£ 



0.05 
0.15 
0.10 
0.12 
0.22 
0.12 
0.025 
0.017 
<0.03 
0.07 
0.05 



(^deduced from the geochemically determined half-life ratio 128 Te/ 130 Te 

( 2 )geochemical result includes all decay modes; other geochemical determinations only marginally 
agree 



(^radiochemical result, again for all decay modes 



Table 2: (3(3(Qv) half-lives in units of 10 26 years corresponding to (m u ) = 50 meV for nuclear matrix 
elements evaluated in the indicated references. 



Nucleus 


Ref.: © 


(80) 


dsij) 


(0) 


(24,0) 


(0) 


48 Ca 


12.7 


35.3 








10.0 


76 Ge 


6.8 


70.8 


56.0 


9.3 


12.8 


14.4 


82 Se 


2.3 


9.6 


22.4 


2.4 


3.2 


6.0 


100 Mo 






4.0 


5.1 


1.2 


15.6 


116 Cd 








1.9 


3.1 


18.8 


iso Te 


0.6 


23.2 


2.8 


2.0 


3.6 


3.4 


136 Xe 




48.4 


13.2 


8.8 


21.2 


7.2 


150 Nd a ) 








0.1 


0.2 




160 Gd a ) 








3.4 







') deformed nucleus; deformation not taken into account 
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Table 3: Best reported limits on Ty 2 - The (m v ) limits and ranges are those deduced by the authors 
and their choices of matrix elements within the cited experimental papers. All are quoted at the 90% 
confidence level except as noted. The range of matrix elements that relate T®Y 2 to (m u ) can be found 
in Table |. 



Isotope Tw 2 (y) (m u ) (eV) Reference 



48 Ca 


> 9.5 x 10 21 (76%) 


< 


8.3 






76 Ge 


> 1.9 x 10 25 


< 


0.35 




© 




> 1.6 x 10 25 


< 


0.33 - 


- 1.35 


© 


82 Se 


> 2.7 x 10 22 (68%) 


< 


5 




© 


100 Mo 


> 5.5 x 10 22 


< 


2.1 






116 Cd 


> 7 x 10 22 


< 


2.6 




© 


128,130fp e 


= (3.52 ± 0.11) xHT 4 
(geochemical) 


< 


1.1 - 


1.5 


© 


128 Te 


> 7.7 x 10 24 


< 


1.1 - 


1.5 


© 


130 Te 


> 1.4 x 10 23 


< 


1.1 - 


2.6 


(ED 


136 Xe 


> 4.4 x 10 23 


< 


1.8 - 


5.2 


(H) 


150 Nd 


> 1.2 x 10 21 


< 


3 




© 



Table 4: The most restrictive {3{3(0v, x) limits. The {g v>x ) limits are those deduced by the authors 
of the cited experimental papers. All are quoted at the 90% confidence level except as noted. The 
total geochemical measured decay rate of 128 Te is used as its f3[3(0v, x) limit. The range of matrix 
elements that relate T^ 2 X to {g u ,x) can be deduced from the entries in Table [| 



Isotope T°^ x (y) 



(fJ, 



Reference 



48 Ca 


> 


7.2 


x 


10 20 


76 Ge 


> 


6.4 


x 


10 22 


82 Se 


> 


2.4 


x 


10 21 


96 Zr 


> 


3.5 


X 


10 20 


100 Mo 


> 


5.4 


X 


10 21 


116 Cd 


> 


3.7 


X 


10 21 


128 Te 


> 


7.7 


X 


10 24 


iao Te 


> 


1.4 


X 


10 21 


136 Xe 


> 


7.2 


X 


10 21 


i50 Nd 


> 


2.8 


X 


10 20 



< 5.3 x 10~ 4 

< 8.1 x 10" 5 

< (2.3 - 4.3) x 10~ 4 

< (2.6-4.9) x 10~ 4 

< 7.3 x 10" 5 

< 1.2 x 10~ 4 

< 3 x 10~ 5 

< (2.6 - 6.7) x 10~ 4 

< (1.3-3.8) x 10~ 4 

< 1 x 10~ 4 



© 



(|100[) 



© 
© 



((TOT]) 
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Table 5: Proposed or suggested future /3f3(0v) experiments, separated into two groups based on 
the magnitude of the proposed isotope mass. The T®J 2 sensitivities are those estimated by the 
collaborators but scaled for 5 years of data taking. These anticipated limits should be used with 
caution since they are based on assumptions about backgrounds for experiments that do not yet 
exist. Since some proposals are more conservative than others in their background estimates, one 
should refrain from using this table to contrast the experiments. The range of matrix elements that 
relate T°Y 2 to (m J) can be found in Table 0. 



Experiment 
COBR Aflllil ) 



Source 
i3U Te 
i50 Nd 

100 Mo 



Detector Description 



Sensitivity to 



TxkF 

2 x 10 25 
4 x 10 24 



DCBA dliir 
NEMO 3^U3|) 



CAMEO ( |114|j~ 
CANDLESmT|) 
CUOREpq) 
EXOp|) 
GEM TO) 
GENIUS^) 
GS0(|11S|,|12C) 



48 Ca 
iso Te 

136 Xe 
76 Ge 
76 Ge 
160 Gd 
76 Ge 
100 Mo 
136 Xe 
136 Xe 



10 kg CdTe semiconductors 

20 kg enr Nd layers between tracking chambers 

10 kg of PP(0u) isotopes (7 kg Mo) with tracking 



> 10^ b 

1 x 10 26 

2 x 10 26 
8 x 10 26 
7 x 10 27 

1 x 10 28 

2 x 10 26 

3 x 10 27 
1 x 10 27 
5 x 10 26 
3 x 10 26 



Majora na(|121 
MOON([l2l)~ 
Xe(p|) 
XMASSQ 



1 t CdW04 crystals in liq. scint. 

several tons of CaF2 crystals in liq. scint. 

750 kg Te02 bolometers 

1 t enr Xe TPC (gas or liquid) 

1 t enr Ge diodes in liq. nitrogen 

1 t 86% enr Ge diodes in liq. nitrogen 

2 t Gd2SiOs:Ce crystal scint. in liq. scint. 
0.5 t 86% segmented enr Ge diodes 

34 t ""'Mo sheets between plastic scint. 
1.56 t of em Xe in liq. scint. 
10 t of liq. Xe 



